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ABSTRACT 

Velocity dispersion measurements of recently discovered Milky Way satellites 
with My > —7 imply they posses high mass-to-light ratios. The expected veloc- 
ity dispersions due to their baryonic mass are ~ 0.2 km s -1 , but values > 3kms _1 
are measured. We perform Monte Carlo simulations of mock radial velocity mea- 
surements of these systems assuming they have mass-to-light ratios similar to 
globular clusters and posses an unidentified binary star population, to determine 
if these stars could boost the velocity dispersion to the observed values. We 
find that this hypothesis is unlikely to produce dispersions much in excess of 
~ 4.5 km s -1 , in agreement with previous work. However, for the systems with 
potentially the smallest velocity dispersions, values consistent with observations 
are produced in 5 — 40% of our simulations for binary fractions in excess of 
fbin(P < lOyrs) ~ 5%. This sample includes the dwarf galaxy candidates that 
lie closest to classical globular clusters in My — space. Considered as a popu- 
lation, it is unlikely that all of these dwarf galaxy candidates have mass-to-light 
ratios typical of globular clusters, but boosting of the observed dispersion by 
binaries from near-zero values cannot be ruled out at high confidence for several 
individual dwarf galaxy candidates. Given the importance of obtaining accurate 
velocity dispersions and dynamical masses for the faintest satellites, it is clearly 
desirable to exclude directly the possible effect of binaries on these systems. This 
requires multi-epoch radial velocity measurements with individual uncertainties 
of <1 kms -1 to identify spectroscopic binaries with orbital velocities of order the 
observed velocity dispersion. 

Subject headings: binaries: general — galaxies: dwarf — galaxies: kinematics and 
dynamics — galaxies: star clusters: general — galaxies: structure — Local Group 
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Introduction 



The internal velocity dispersions of Milky Way (M W) dwarf galaxies with My < —8 



lie in the range 6 < a b s < llkms 1 ( lAaronson 



1983). These values are large given the 



observed luminosities and surface brightnesses, and prompted early research to consider 



unidentified spectroscopic binary stars as a 



velocity dispersions (e.g., 



Mateo et al. 



1993 



potential explanation for the unexpected 



Olszewski et al. 



1996 



Hargreaves et al. 



y hig h 



1996|). 



If many of the targeted stars are binaries, then the radial velocity that is measured is a 
superposition of their intrinsic radial velocity within the galaxy and their binary orbital 
velocity. Simulations of plausible binary populations showed that, while some inflation of 
the dispersion may be expected due to binaries, this effect is generally insufficient to explain 



the high ve 



Minor et al. 



o city d ispersions that are measured for the dwarf galaxies. More recently, 



(120101 ) develop a methodology to determine the contribution to the velocity 
dispersion of a dwarf galaxy from a binary population with a given set of parameters, for a 
well sampled dataset. 

Here we revisit the role of binary stars in boosting the observed velocity dispersion of 
intrinsically low-mass systems. We are motivated by the discoveries of MW satellites with 
considerably lower luminosity (My > —7) than previously known dwarf galaxies, a regime 
that was not considered in the original papers on this subject. Some of these systems have 
few bright stars suitable for spectroscopic follow-up, and a few have measured velocity 
dispersions as low as a obs ~ 3kms _1 (implying mass-to-light ratios T v = M/L v > 100). 
In Section 2, we summarize the relevant properties of these satellites, and we discuss how 
we generate realistic binary populations. In Section 3, we use Monte-Carlo simulations to 
investigate the effect of binary stars on velocity dispersion measurements of these systems 
assuming they have a near-zero intrinsic dispersion. 
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2. Preliminaries 
2.1. The faint population of Milky Way satellites 

Fourteen new MW sate 



2006 



2006 



Zucker et a 



Irwin et al 



2006a 



2007 



lites discovered since 2004 



Belokurov et al 



Walsh et al 



2006 



2007 



Willman et al. 


2005: 


Willman et al. 


, 2008 




2009; 


Sakamoto & Haseeawa 



20071 ) have been shown to possess velocity dispersions 



significantly higher than expe cted if they have mass-to-light ratios typical of globular 



clusters in virial 


equilibrium 


L f 


Martin et al. 


2007 


Koch et al. 


2009 


• 


Walker et al. 


2009 




Aden et al. 



Simon & Geha 


2007 




Geha et al. 


2009 


2009 


> 


Belokurov et al. 


2009 


) . Their 



luminosities span more than five magnitudes, with the brightest systems (Canes Venatici I 
and Leo T) of comparable luminosity to Draco (the faintest MW dwarf galaxy known 
previously) . 

Table 1 summarizes the relevant structural and kinematic data for recently discovered 
faint MW satellites. We define this as My > —7, and note that this sample includes three 
systems (Ursa Major I, II and Bootes I) with velocity dispersions of the same order as those 
of previously known dwarfs. In addition to the observed velocity dispersion (a b s ), we list 
the number of putative member stars each measurement is based on (n stars ), and the median 
velocity uncertainty (v err ) of those members. Leo V has two dispersion measurements, a 
lower central value using five stars and a larger global value using a n additional two sta rs 



located at larger radius. Here, we quote the larger value. Further, 



Walker et al. 



(|2009|) 



do not rule out that Leo V may be a star cluster since its dispersion is 



reso 
and 



ved. We do not in clude Bootes III (cr (, s ~ 14kms 1 ) in Table 1; iGrillmairl ( 120091 ) 



presently 



jarely 



Carlin et al. 



(120091 ) both suggest that Bootes III is a dwarf galaxy undergoing tidal 



^Belokurov et al.l ( 120101 ) recently discovered Pisces II, a candidate MW dwarf galaxy that 
does not yet have a velocity dispersion measurement. 
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disruption. While we confine our analysis to the MW satellite population, we note that 
a similar analysis can also be applied to other Local Group satellites, suc h as the M31 
ite system, for whi ch kinematic data is now becoming available (e.g 



satel 



2010 



Kalirai et al. 



Collins et al. 



20101 ). 



For each satellite in Table 1, we calculate what their properties would be if they possess 
structural and kinematic characteristics typical for globular clusters (whose mass-to-light 
ratios are consistent with expectations from population synthesis models, without dark 
m atter) . Specifically, we derive their concentration, c, and core radius, r c using the relations 



of 



McLaughlin! (120001 ) . and derive the corresponding intrinsic velocity dispersion, <jj nt , using 



0~int 



27rGapM\* 



9ur r 



(1) 



assuming M = TyLy where Ty = 1.45 ± 0.1 in solar units ( McLaughlin 



constants b>, a and p ar e calculated by interpolation of the values given in 



2000) 



King 



. The 



(119661 ) and 



Peterson fc King! ( 119751 ). Results are shown in Table 2. The uncertainty on Oi n t corresponds 
to the uncertainty on Ty. The derived values for a int are ~ 0.2kms _1 , generally at least 
an order of magnitude smaller than a b s and v err . 



2.2. Calculating the velocity dispersion 

We will derive velocity dispersion estimates from mock radial velocity catalogs for 
these satellites by attempting to reproduce the observational techniques used to measure 
(fobs. In general, this involves identifying stars belonging to the system and measuring their 
mean velocity and dispersion. Here, we know a priori that our stars are members, but we 
will not necessarily use all the stars in our mock catalog to measure the dispersion. For 
example, some short-period binaries may have very different velocities from the mean of the 
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sample, and they may not be observationally identified as members of the satellite. Instead, 
we sigma-clip our samples at 3a in order to mimic observations as closely as possible (e.g., 



Simon fc Geha 



2007 



Martin et al. 



20071 ). At most, this removes only 1 or 2 stars from 
the final sample. All of the putative member stars identified in the actual observations 
summarized in Table 1 lie within 3a of the mean velocity, and so this seems reasonable. 

Several of the radial velocity studies listed in Table 1 calculate the velocity dispersion 
(op) and mean velocity (u) of the putative members by maximizing the natural logarithm 
of the probability function, 



ln(p) 



1 (v 

_Y^ n atar V 

2 



u) 2 



1 « 



+ a 2 ) 

err 1 pj 



T^stars 



ln(27r) 



(2) 



( IWalker et al. 



20061 1. where Vi >err is the measurement uncertainty on the measured velocity, 
Vi, of the i th star. We use this technique to calculate the velocity dispersion of our mock 
datasets for every satellite. We test this method on the observed datasets, and derive the 
same cen tral values as the publi shed results, with the exception of Segue II a nd Leo V. For 



Segue II 



Belokurov et al 



Walker et al. 



(120091 ) use a different likelihood analysis; for Leo V 
( 120091 ) use a similar technique to Equation 2 but incorporate the probability that a star is 
a member of the system. To ensure consistency, we recalculate the velocity dispersions for 
Segue II and Leo V using Equation 2 and use these new values in our subsequent analysis. 
These are inc l uded in Table 1 in parentheses, with the uncertainties calculated following 



Walker et al 



( 120061 ) 
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2.3. The binary star population 

The observed radial velocity of the primary companion in a binary star system is given 

by 

2ixai sin i 

V r — [cos(B + oj) + ecosw] . (3) 

P(l — e 2 )2 

ai corresponds to the semi-major axis of the relative orbit of the primary; P is the period 
of the orbit; e is the eccentricity; i is the orbital inclination to our line-of-sight; 9 is the 
phase of the binary measured from periastron; oj is the longitude of periastron; mi and m 2 
are the masses of the primary and secondary, respectively, m = m 1 + m 2 and the mass 
ratio q = 777.2/ 'vn,\. For an equal- mass binary system on a circular orbit with m = 1 M Q and 
P = 1 yr, the maximum radial velocity that can be observed is V r ^ max ~ ±15kms -1 . For 
P = [10, 100, 1000] yrs, V r , max ~ [±7, 3, 1.5] kms" 1 . 

We now summarize the adopted distributions for each of the variables in Equation 3. 
At the present time, the relevant physical parameters are entirely unknown for the stellar 
systems under discussion, and so we must be guided in our choice by the findings from other 
stellar populations, particularly the solar neighborhood (SN) and halo globular clusters. 

2.3.1. Mass distributions 

We set the mass of the primary to be m\ = 0.8 M ; this is the approximate minimum 
mass of a star that has had time to evolve on to the giant branch. The studies summarised 
in Table 1 target either giants or stars near the main-sequence turn-off, and so this seems 
reasonable. 
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Duquennoy fc Mayorl (Il99ll ) show that the mass ratio distribution for binary stars in 
the SN is well approximated as log-normal, 



dN 
dq 



oc exp 



-\2 



(q-q) 

2a\ 



(4) 



q = 0.23 and a q = 0.42. We set q min = 0.1, so that the minimum mass of the secondary 
approximately corresponds to the limiting mass for hydrogen burning. 



2.3.2. Period distribution 



We use two different binary period distributions. We first follow 



Duquennoy fc Mayor 



(I199ll ). who show that the distribution for SN binaries is described by a log- normal, 



dN 



dlog 10 P 



oc exp 



(log w P - logipP)' 
In 2 

Z<J log 10 P 



(5) 



logioP = 4.8 and o~i g 10 p = 2.3, where P is measured in days. We also explore a uniform 
distribution in log w P. 

The maximum possible period of a binary star in a stellar cluster, P ma x, is set physically 
by the maximum separation b eyond which the binary becomes unbound due to stellar 

(jl984J) derive the semi-major axis at which this occurs, 



interactions within the cluster. 



Hills 



12 AAU 



mi + 1^2 



10 kms 



-i 



a 



(6) 



where a corresponds to the intrinsic velocity dispersion of the cluster. For the values of 
cr iri t in Table 2, this corresponds to P max ~ 10'yrs, or V r . max ~ 0.14 kms ^, of order o~i n t. 
However, such systems, if present, are likely very rare (following Equation 5) and cannot 
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contribute to inflating the velocity dispersion by more than a small fraction of a kms 1 . 

Constraints on P m i n can be derived physically for red giant branch (RGB) stars from 
mass transfer considerations, 



h(q) 



0-38 < h(q) < 0.59 



(7) 



where R\ is the radius of a star on the RGB and h(q) is a slowly varying function that 
ume radius of the Roche lobe in terms of the orbital separation (see §V of 



gives the vo 


ume 


Prvor et al. 


1988) 



and so we conservatively adopt (a min ) ~ 1 AU. Thus, P m in ~ 1 y r typically, although values 
as low as 0.1 yr are plausible, particularly for stars fainter than the tip of the RGB. 



2.3.3. Eccentricity distribution 



Two eccentricity distributions are used. Following 
thermal distribution in eccentricity, 



Heggid ( 119751 ). we first assume a 



dN 
de 



oc 2e . 



In the second instance, we assume circular orbits throughout (e = 0) 



2.3.4. Angles 

The probability of observing a binary system at an inclination i is proportional to 
sini. The probability of observing a binary system with a given orientation of its major 
axis to the line of sight, u, is uniform between and tt. The probability of observing a 
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binary system at a given phase, 6, is proportional to the inverse of the angular velocity, 
^ _1 (m,P, e). Thus, for circular orbits, the probability of observing a binary star at a given 
6 is uniformly distributed. For elliptical orbits, the binary star is more likely to be observed 
near apastron. 



2.3.5. Binary fractions 



The fraction of stars that are binaries in the systems under discussion (fun) is 
presently entirely unconstrained, and there is some uncertainty as to the binary fraction 



i n oth er stellar populations. For example, in the solar neighborhood, 



Duquennoy fc Mayor 



( 1l99ll ) find a n overall b i nary fraction of ~2/3 for stars near the middle of the main 



sequence and 



Hut et al 



(119921 ) conclude that the binary fraction in globular clusters is only 



slightly smaller t han t hat in the Galactic disk. Few constraints exist in dwarf galaxies; 



Olszewski et al. 



(jl996[ ) estimate a binary fraction in Draco and Ursa Minor of 0.2 — 0.3 per 



decade in period, for periods around 1 yr, that is potentially 3 — 5 times 



solar neighborhood, but which is also highly uncertain. However, 



ligher than in the 



Ladal ( 120061 ) argue that 



most stars are in fact single, with two-thirds of main sequence sy stems in the Gala ctic disk 



Cote et al. 



fjl996T ) have 



composed of single stars. In a small sample of globular clusters, 
estimated a binary fraction of ~ 0.01 — 0.0 8 per decade of period (with significant error 



bars) , and other studies such as 



Davis et al 



(120081 ) in NGC 6397 imply a present-day binary 



fraction of order a perc ent. Formation mo dels of glo bular clusters can c reate primordial 



binary fractions of 5% (IHurley et al. 



20071 ) to 100% (llvanova et al.ll2005h . In short, it is 



presently unclear what the primordial or current binary fraction of the systems listed in 
Table 1 will be. 

Given the lack of constraints on /^ n , we examine three different possible normalizations, 
fbin{P — lOyrs) = 0.05, 0.15 and 0.30. The latter normalization is high, but is similar in 
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magnitude to that estimated by 



Olszewski et al. 



(Il996[ ). The period distribution must be 



truncated at some upper limit else overall binary fractions > 1 can be obtained depending 
on the shape of the assumed distribution; since we are only interested in binaries numerous 
enough and with short-enough periods to contribute to <r b s , we do not consider periods in 
excess of 10 — 100 yrs. 



3. Analysis 
3.1. Method 

We simulate four different binary populations with different assumptions about the 
period and eccentricity distributions: (i) log-normal in P, thermal in e; (ii) log-normal in 
P, e = 0; (iii) uniform in logioP, thermal in e; (iv) uniform in logioP, e = 0. In each case, 
we generate 10 6 binary star systems and output the results to a file. 

For each of the binary populations considered, we generate mock velocity datasets for 
each satellite using the following procedure: 

1. For i = 1, ...n s tars, we select an intrinsic velocity Vi for each star from a Gaussian 
distribution with width a int ; 

2. For each star, we randomly designate it as single or binary according to fun\ 

3. For stars that are designated binary, we randomly assign them one of the 10 6 binary 
orbital properties generated previously (ensuring P m i n < P < P ma x)- A new value for 
v i is obtained by summing the previous value and the orbital radial velocity of the 
primary; 

4. For each star, we assign it an observational uncertainty from the list of n star 
observational uncertainties (v err ^) that were actually obtained through observations 
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for that satellite. These are available in the corresponding papers referenced in Table 2 
or were provided by the authors (J. Simon & M. Geha, private communication). Vt is 
then modified by an error-term selected from a Gaussian distribution of width v err /, 

5. Once n star mock radial velocities are generated according to (1) — (4), we make a 
first estimate of the mean (— * -) and standard deviation (\ — * . ). All stars 

v n s tar ' v V nstar — 1 ' 

that are more than 3 standard deviations from the mean velocity are not considered 
further; 

6. For the sigma-clipped sample of stars, we calculate and record the final values for the 
mean velocity and velocity dispersion using Equation 2; 

7. Steps (1) — (6) are repeated 10 000 times for each satellite. 

We record the percentage of times the final velocity dispersion estimate is greater than 
or equal to the lower la error bound of a obs . These percentages are listed in Table 3 for 
each of the binary populations considered, for different values of P m i n and P max , assuming 
different binary fractions. Figure 1 shows examples of the distribution of measured 
velocity dispersions from these simulations for Ursa Major I (left panel), Segue (middle 
panel) and Leo IV (right panel). In each panel, we show the results corresponding to 
fbin(P < lOyrs) = 0.15, a log-normal binary period distribution (1 < P < 100 yrs) and 
circular orbits (column 7 of Table 3). 

3.2. Results 

Table 3 shows that the binary populations considered here are very likely insufficient to 
explain the velocity dispersion measurements for Ursa Major I, II and Bootes I (the three 
systems with largest a b s ; see the left panel of Figure 1), and similarly for Coma Berenices 
and Canes Venatici II (unless there is a high binary fraction in these systems). 
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For the remaining systems in Table 3 (Segue I, II, Willman I, Bootes II, Leo IV, V and 
Hercules), there is a modest, but non- negligible, chance (> 5%) that the measured values 
of a b s may be due to binary stars, unless the binary fraction is < 5%. For higher binary 
fractions, the probability that binary stars can produce a Q b s can increase to > 20%. All of 
the satellites most affected have a Q b s < 4.5km s _1 (with the exception of Bootes II, but here 
the large observational uncertainty allows a value as low as 3.1km s" 1 ). The middle panel 
of Figure 1 shows an example velocity dispersion distribution for Segue, a fairly typical 
case. Leo IV is the satellite potentially most affected by binaries (right panel of Figure 1), 
and inspection of Table 1 shows that this is a result of a relatively low value of a D b s in 
combination with a relatively large v err . 

Given that we do not know the actual binary fractions in any of the systems listed in 
Table 1, or the distribution of their properties, we conclude that the present observational 
data on the satellites with lowest a Q b s are unable to exclude at high confidence that their 
velocity dispersions may result from binary stars in intrinsically low-mass systems. On the 
other hand, a b s for Ursa Major I is extremely unlikely to be due to binary stars, and results 
for Ursa Major II, Bootes I, Coma Berenices and Canes Venatici II are unlikely to be able to 
be explained by binary stars unless these systems possess high binary fractions. Considered 
as a population, it is unlikely that all these systems have mass-to-light ratios typical of 
globular clusters, but boosting of the observed dispersion by binaries from near-zero values 
cannot be ruled out at high confidence for several individual dwarf galaxy candidates. 



I n line with earlier stud i es that dealt with higher velocity dispersion dwarfs (IMateo et al. 



1993 



Olszewski et al. 



1996 



Hargreaves et al. 



19961 ). we find binary stars cannot boost 



velocity dispersions to values much in excess of 4.5km s 1 . Our results are also consistent 

i — — i 

with those of Minor et al.l (120101 ). who examine the effect of binaries on dwarf galaxies with 
G%nt ^ 4kms _1 , which is an order of magnitude larger than the values of Oi nt considered 
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here. 

The binary stars that contribute to boosting the velocity dispersion from a int to a obs 
must have current orbital radial velocities of order a Q b s - In our simulated data, we find that 
the average yearly change in the orbital radial velocities of such stars is ~ 1 — 3kms _1 . To 
identify such stars, multi-epoch measurements (with a cadence of a year) with individual 
velocity uncertainties < 1km s -1 would be required; these observations will then directly 
determine the possible influence of binary stars on cr b s . 

4. Summary 

We have examined the possibility that unidentified binary stars may inflate the 
velocity dispersions measured for some low-mass MW satellites. If such objects were to 
have mass-to-light ratios similar to globular clusters, then they would be expected to have 
intrinsic velocity dispersions of 0.1 < a int < 0.3km s" 1 — significantly below the precision 
of existing radial velocity measurements in these systems. Our Monte Carlo simulations 
reveal that, although binaries are unlikely to explain satellite velocity dispersions much in 
excess of ~4.5 km/s, boosting of the observed dispersion by binaries cannot be ruled out 
with high confidence for some faint satellite candidates. Figure 2 — which plots r\ versus 
My for MW globular clusters, dwarf galaxies and dwarf galaxy candidates — shows that 
it is those dwarf candidates with r h < 70 pc (i.e. lying nearest to the globular cluster 
population) that are the most susceptible to contamination by spectroscopic binaries. 
Given the importance of obtaining accurate velocity dispersions and dynamical masses 
for the faintest MW satellites, it is clearly desirable to exclude directly the possible effect 
of binaries on these systems by obtaining multi-epoch radial velocity measurements with 
individual uncertainties of <lkms~ 1 . 
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Table 1: Observed parameters and relevant kinematic data for faint MW satellites. 
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Table 2: Derived characteristics of faint MW satellites assuming Y v = 1.45. 
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Table 3: Simulation results for each satellite assuming fbi n (P < 10 yrs) = 5%, 15% and 
30% (upper, middle and lower sections, respectively), for different period and eccentricity 



assumptions. Numbers refer to the percentage of simulations that generate a velocity dis- 
persion that is greater than or equal to the la lower error-bound on the observed velocity 
dispersion. 
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Fig. 1. — The simulated velocity dispersion distribution for Ursa Major I, Segue, and Leo IV, 
assuming fbi n (P < lOyrs) = 15%, a log-normal period distribution with 1 < P < 100 yrs and 
circular orbits (column 7 of Table 3). The dotted line shows the assumed intrinsic velocity 
dispersion, and the dashed lines show the observed velocity dispersion with la error bounds. 
Ursa Major I is incompatible with an intrinsically low mass system, but this hypothesis 
cannot be excluded at high confidence for Segue and Leo IV. 
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dwarfs whose velocity dispersions are most susceptible to boosting by spectroscopic binaries 
are shown as filled squares. These include all dwarf candidates with < 70 pc that lie 



closest to the globular cluster population. 



